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1 Introduction 



A dispersionless limit of PDEs is taken by rescaling the independent variables X a — ► X a /e, and 
taking the limit e — > 0. This is a delicate procedure, as the limit of the solutions of a given PDE does 
not usually correspond to solutions of the limiting dispersionless equation. Moreover, inequivalent 
PDEs may have the same dispersionless limit, so the problem 

• Recover the original PDE from its dispersionless limit 

is, of course, ill posed. Some progress can nevertheless be made if the dispersionless equation is 
integrable, and one insists that its dispersive analogue is also integrable. In the next section we 
shall explain how dispersionless limits of solitonic PDEs are equivalent to the WKB quasi-classical 
approximation of the associated linear problems. This suggests that the reconstruction of the dispersive 
solitonic system should involve a quantisation of some kind. 

Such a quantisation procedure has been developed in the seminal work of Kupershmidt [15]. This 
procedure is based on the Moyal product, and works well if the Lie algebra underlying the dispersionless 
linear problem is the algebra sdiff(S 2 ) of divergence-free vector fields on a 2-surface S 2 . This is the 
case for the dKP and SU(oo) Toda equations in 2 + 1 dimensions. Similar progress can also be made 
in higher dimensions and, indeed, one of us has constructed integrable deformations of Plebanski's 
first heavenly equation [29] by replacing the underlying Poisson bracket with the Moyal bracket. 

The idea of deforming integrable systems while retaining the integrability of the resulting equation 
has now been studied from a number of different points of view: 

• Takasaki studied properties of the deformed heavenly equations and described how solutions may 
be described in terms of a Riemann-Hilbert splitting in a Moyal algebra valued loop group [33]. 
Extensions of this led to Moyal-KP hierarchies [34] and deformations of the self-dual Yang- 
Mills equations [35]. The deformed Riemann-Hilbert procedure was recently fully developed by 
Formanski and Przanowski [9, 10]. 

• Nekrasov and Schwarz introduced instantons on noncommutative space-time [22]. This led to 
the development of noncommutative soliton equations. These may be viewed as a deformation 
of the standard, commutative, soliton equations. Many of these may be studied as reductions of 
the noncommutative self-dual Yang-Mills equations [12, 16]. 

• Associated to any Frobenius manifold is a hierarchy of integrable equations of hydrodynamic 
type. Integrable deformations of these equations arise naturally when one studies the genus 
expansion in the corresponding topological quantum field theories [3]. 

In the present paper deformations of multidimensional integrable systems are based on the algebra 1 
diff(S), the Lie algebra of vector fields on S, where S = 5 1 or R. It turns out, however, that 
this algebra admits no non-trivial deformations [17]. However an alternative method of deforming 
these integrable systems may be developed. This method is based on the approach of Ovsienko 
and Rogers [24] where a homomorphism from diff(S) to the Poisson algebra on T*E can be used to 
construct non-trivial deformations. We shall use this idea to construct integrable deformations of 
various equations associated to the algebra diff(S) . 

1 In the remainder of this paper the superscript, denoting the dimension of the manifold, will be dropped. 
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It should be pointed out that the original, undeformed, equations have a natural interpretation 
in terms of twistor theory, via the non-linear graviton construction and its variants. It would seem 
desirable to develop a 'deformed' version of twistor theory that would encode solutions of the deformed 
equations as some sort of deformed holomorphic conditions. This idea was what was behind the 
paper [31], but the problem remains open (though see [20] for some ideas on how Nijenhuis structures 
may be deformed). 



2 Dispersionless limit in 2 + 1 dimensions 

Certain dispersionless integrable systems can arise from solitonic systems in a following way: Let 

Q Qn Qn-1 

a (bx) = 555 + «'(*->55i=T + - + «.<*•>. 

g Qm Qm-l 

be differential operators on R with coefficients depending on local coordinates X a = (X, Y, T) on R 3 . 
The overdetermined linear system 

admits a solution %f(X,Y,T) on a neighbourhood of initial point (X,Yq,Tq) for arbitrary initial data 
^(X, Yq,Tq) = f(X) if and only if the integrability conditions fy-r = ^ty, or 

A T -B Y + [A,B] =0 (2.1) 

are satisfied. The nonlinear system (2.1) for oi, . . . , a n , b±, . . . , b m can be solved by the inverse scat- 
tering transform (1ST). Integrable systems which admit a Lax representation (2.1) will be referred to 
as solitonic, or dispersive. 

The dispersionless limit [37] is obtained by substituting 

e-£-, V(X a )=exp(iP(x a /e)) 1 



Q X a Q x a ' 

and taking the limit e — ► 0. In the limit the commutators of differential operators are replaced by 
the Poisson brackets of their symbols according to the relation 

d k , , , r . „, dAdB dAdB 



¥ — (TM**, IA,B]^——-—— = {A,B}, A = ^, 



dX k ' ' d\ dx dx d\ 

where A, B are polynomials in A, with coefficients depending on x a = (x, y, t). The dispersionless limit 
of the system (2.1) is 

A t - B y + {A, B} = 0. (2.2) 

Nonlinear differential equations of the form (2.2) are called dispersionless integrable systems. One 
motivation for studying of dispersionless integrable systems is their role in constructing partition 
functions in topological field theories [14]. 
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A natural approach to solving (2.2) would be an attempt to take a quasi-classical limit of the 1ST 
which linearises (2.1). This does not yield the expected result, as the quasi-classical limit of the Lax 
representation for (2.1) is the system of Hamilton-Jacobi equations 

ip y = A{ip x ,x a ), Vt = B(i(> x ,x a ), 

with 'two times' t and y, and the initial value problem for (2.2) would require a reconstruction of a 
potential from the asymptotic form of the Hamiltonians. This classical inverse scattering problem is 
so far open. 

There are alternative methods of solving (2.2) [13, 8, 32, 6]. In particular the mini-twistor approach 
of [6] works as follows: The system (2.2) is equivalent to the integrability [L, M] = of a two- 
dimensional distribution of vector fields 

L = d t - B x d x + B x d x , M = d y - A\d x + A x d\ (2.3) 

on R 3 x MP 1 . Assume that L,M are real analytic, and complexify M 3 to C 3 . The mini-twistor space 
Z is the two complex dimensional quotient manifold 

Z = C 3 x CP7(£,M), AeCP\x a eC 3 . 

That is to say that the local coordinates on Z lift to functions on C 3 x CP 1 constant along L, M. 

The mini-twistor space is equipped with a three parameter family of certain rational curves. All 
solutions to (2.2) can in principle be reconstructed from a complex structure of the mini-twistor space. 

In fact the twistor approach outlined above is capable of solving a wider class of equations. We 
shall therefore generalise the notion of the dispersionless integrable systems by allowing distributions 
of vector fields more general than (2.3). The derivatives A\,A X ,B\,B X of the symbols (A,B) of 
operators can be replaced by independent polynomials A± , A 2 , B\ , B2 in A with coefficients depending 
on (x, t/, t) 

L = 8 t - B x d x + B 2 d x , M = d y - A 1 d x + A 2 8 x . (2.4) 

If A\,B\ are linear in A and A 2 , B 2 are at most cubic in A then the rational curves in Z have normal 
bundle 0(2) (the line bundle over CP 1 with transition functions A~ 2 from the set A / 00 to A / i.e. 
Chern class 2) and the three-dimensional moduli space of such curves in Z can be parametrised by 
(x,y,t). Allowing polynomials of higher degrees would lead to hierarchies of dispersionless equations. 
We take the integrability of this generalised distribution (2.4) as our definition of the dispersionless 
integrable system. The definition is intrinsic in a sense that it does not refer to an underlying dispersive 
equation. 

3 Diff(5' 1 ) dispersionless integrable systems 

In this section two integrable systems associated with the gauge group Diff(5' 1 ) will be given. The 
first has been extensively studied in [26, 8, 5, 19, 7, 18, 25], so only a new gauge theoretic description 
will be given - the reader is referred to these earlier papers for more details. The second system, which 
arises from a Nahm-type system, is new and this system is discussed in more detail. 
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3.1 A (2 + 1) dimensional dispersionless integrable system 

An example of a dispersionless system which is integrable in the sense of the outlined twistor corre- 
spondence is given by the following distribution 

L = dt — wd x — \d y , M = d y + ud x — Xd x . (3.1) 

A linear combination of this distribution leads to a special case of (2.4) with A2 = B2 = 0. Its 
integrability leads to the pair of quasi-linear PDEs 

fH + w y + uw x — wu x = 0, u y + w x = 0, (3-2) 

for two real functions u = u(x, y,t),w = w(x, y, t). This system of equations has recently been studied 
in [26, 8, 5, 19, 7, 18, 25] in connection with Einstein-Weyl geometry, hydrodynamic chains and 
symmetry reductions of anti-self-dual Yang-Mills equations. From the twistor point of view (3.2) is 
invariantly characterised [5] by requiring that the mini-twistor space Z fibres holomorphically over 
CP 1 . The second equation can be used to introduce a potential H such that u = H x ,w = —H y . The 
first equation then gives 

Hxt — H yy + H y H xx — H x H xy = 0. (3-3) 

The system (3.2) arises as a symmetry reduction of the anti-self-dual Yang Mills equations in 
signature (2, 2) with the infinite-dimensional gauge group Diff(S) and two commuting translational 
symmetries exactly one of which is null [7]. This combined with the embedding of SU(1, 1) C Diff(S) 
gives rise to explicit solutions to (3.2) in terms of solutions to the nonlinear Schrodinger equation, and 
the Korteweg de Vries equation [7]. 

The Lie algebra of the group of diffeomorphisms Diff(S), where £ = S 1 or R, is isomorphic to the 
infinite-dimensional Lie algebra of functions on S with the Wronskian 

< f,9>-= f9x~ fx9 (3-4) 

as the Lie bracket, where f,g G C°°(S), and x is a local coordinate on S. An alternative gauge- 
theoretic interpretation can be given to (3.2): Observe that the first equation in (3.2) can be interpreted 
as the flatness of a gauge connection on R 2 , where the gauge group is Diff(S). Indeed, choose local 
coordinates (t,y) on R 2 and consider A G A X (R 2 ) <g> C°°(S) of the general form 

A = —wdt + udy, 

where u, w : R 2 — ► C°°(E) depend on (x,y,t). 
The flatness of this connection yields 

dA+ AAA = (u t + w y + <u,w >)dt A dy = 0, 

as claimed. Therefore the connection is a pure gauge, and can be written as A = g~ 1 dg, where 
g = g(x, y, t) G Map(R 2 , Diff (£)), and 

w = -g~ x gu u = 9~ X 9y 



5 



The second equation in (3.2) yields the following system 

i^ 1 9y)y - {g^Qtjx = 0, (3.5) 
where g = exp (A) is a finite diffeomorphism of S, and terms like g~ x gt should be understood as 

g~ l g t = A t - <A,A t >+^<A,<A,A t »+... . 
3.2 A (3 + 1) dimensional dispersionless integrable system 

In this section we shall present another example of an integrable system associated to the Lie algebra 
of Diff(S' 1 ). We shall first write it as a Nahm system 

e» = -e ijk [ej, e k ], i = 1, 2, 3 (3.6) 

where are vector fields on an open set in R 4 given by 

and (x,y J ) are local coordinates. Rewrite (3.6) as 

d x Ni + e ijk djN k - ^e ijk < Nj , N k >= 0. (3.7) 

We shall now discuss the origin and possible applications of (3.7) 

1. Any solution to (3.7) defines a hyper Hermitian conformal structure represented by the metric 

g = n 2 + 6 ij dy i dyi, (3.8) 

where 

n = dx + Nidy\ 

The three complex structures = 1,2,3 satisfying the algebra of quaternions 
are given by 

Ii(n) = dy\ 

These formulae together with the algebraic relations satisfied by Ij determine the complex struc- 
tures uniquelly, e.g. 

Ii(dy J ) = -<%(n) + e ijk dy k . 

One way to impose integrability of the complex structures is to use the explicit form of the 
complex structures on the basis (dy 1 , dy 2 , dy 3 , n) and demand that the space A^ 1,0 ) is closed 
under exterior differentiation. We begin by defining a basis of self-dual 2-forms 

S' = nA dy* + ^e ijk dy j A dy fe . 
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The integrability of the complex structures is then equivalent to the anti-self-duality of the 
two-form dn: 

S i Adn = 0. (3.9) 

This condition is equivalent to (3.7). 

A dual formulation leads to the Lax pair of vector fields, which is a special form of the hyper- 
Hermitian Lax pair [4, 11, 2]. To see it set = d x and define complex vector fields 

w = ei — ze2, z = e3 — ie^. 

The system (3.7) is equivalent to the commutativity of the Lax pair 

[w-Az,z + Aw] =0 (3.10) 

for all values of the parameter A. 

2. The system (3.7) arises as a symmetry reduction of the anti-self-dual Yang Mills equations on 
R 4 with the infinite-dimensional gauge group Diff(S' 1 ) and one translational symmetry. In fact 
any such symmetry reduction is gauge equivalent to (3.7). To see it consider the flat metric on 
R 4 which in double null coordinates w = y 1 + iy 2 , z = y 3 + iy A takes the form 

ds 2 = dzdz + dwdw, 

and choose the volume element dw A dw A dz A dz. Let A £ T*R ®g be a connection one- 
form, and let F be its curvature two-form. Here q is the Lie algebra of some (possibly infinite 
dimensional) gauge group G. In a local trivialisation A = A^dy^ 1 and F = (l/2)F flu dy tl A dy u , 
where F^ v = D v ] takes its values in q. Here = + A^ is the covariant derivative. The 
connection is defined up to gauge transformations A — ► b~ l Ab — b~ 1 db, where b € Map(IR , G). 
The ASDYM equations on A^ are F = — * F, or 

F — D F — \- F - = F — = f) 

1 wz — u ) - 1 WW ~ - 1 ZZ u ) A wz u - 

These equations are equivalent to the commutativity of the Lax pair 

L = D W -XD-, M = D z + XD W (3.11) 
for every value of the parameter A. 

We shall require that the connection possesses a symmetry which in our coordinates is given by 
d/dy 4 . Choose a gauge such that the Higgs field A^ is a constant in q. Now choose G = Diff(S' 1 ), 
so that the components of the one form A become vector fields on S 1 . We can choose a local 
coordinate x on S 1 such that A^ = d x , and Ai = —Nid x , where Ni = Ni(x,y^) are smooth 
functions on R 4 . The Lax pair (3.11) is identical to (3.10) and the ASDYM equations reduce to 
the first order PDEs (3.7). 

3. Example. An ansatz N(x,y 3 ') = / (x)A(y- 7 ), where N = (Aq, N2, Ar 3 ) T , reduces (3.7) to a pair of 
linear equations 

f = cf, cA + V A A = 0, 
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where c is a constant. If c = then N may be absorbed into a redefinition of the coordinate 
x in the metric (3.8). Therefore we assume c / 0. We set c = 1 by rescaling y 3 and solve for 
/ = exp (x) reabsorbing another constant of integration into A. Now define a new coordinate 
p = exp (—x). Rescaling the metric (3.8) yields 

g = pdy 2 + p- l {dp- A-dy) 2 . (3.12) 

This metric is hyper-hermitian iff the vector A(y J ) satisfies the Beltrami equation 

A + VAA = 0. (3.13) 

This is a slight improvement of the result of [36] where it is claimed that (3.12) is ASD iff (3.13) 
holds. 

The Beltrami equation implies that A is divergence-free and satisfies AA + A = 0, where 
A = V 2 is the scalar Laplacian on M 3 acting on components of A. Existence of solutions of 
equation (3.13), at least in the analytic case, can be proved by an application of the Cartan- 
Kahler theorem (c.f. Example 3.7 in Chapter III of [1]). 

4. The system (3.7) can be put in the hydrodynamic form 

a r N = McurlN, 

where 

/ 1 N 3 -N 2 \ _1 
M = - -N 3 1 JVi . 
V N 2 —Ni 1 / 

A different analytic continuation of (3.7) can be obtained at the level of the hyperHermitian 

geometry. This comes down to looking for conformal structures (3.8) of signature + H . To 

achieve this we regard yi,y 2 ,N±, N 2 as imaginary, and define 

Y 1 = iy 1 , Y 2 = iy 2 , Y 3 = y 3 , Mi = iNi , M 2 = iN 2 , M 3 = N 3 . 

The desired system for Mi = Mi(Y J , x) arises from (3.7). 

Clearly there are many further properties of these dispersionless systems that may be studied. We 
now turn our attention to the construction of non-trivial integrable deformations of equations (3.2) 
and (3.7). 

4 Dispersive deformations 

Given a dispersionless integrable system it is natural to ask whether it arises as a limit of some 
dispersive (or solitonic) system. One would expect the reconstruction of a solitonic system to involve 
a quantisation of some kind, because taking a dispersionless limit of (2.1) was equivalent to a quasi- 
classical limit of the wave function ^(X a ). This is indeed the case, and the paradigm example is 
provided by the connection between the Kadomtsev-Petviashvili (KP) equation and its dispersionless 
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limit dKP. One can reconstruct KP from dKP by expressing the latter in the form (2.2), and replacing 
the Poisson brackets by the Moyal bracket [15, 32, 30, 31, 27]. The infinite series involved in a Moyal 
product truncates in this case, because the symbols A and B are polynomials in momentum A. The 
deformation parameter can then be set to one, and removed from the construction. It seems however 
that this beautiful example is rather exceptional, and that the reconstruction of dispersive systems (if 
at all possible) is in general non-unique and can lead to systems which involve a formal power series. 

Generalising the definition of the dispersionless systems to nom-Hamiltonian distributions (L, M) 
like (3.1) makes things even worse, as the Poisson bracket is not present, and the connection with 
known quantisation procedures of a classical phase-space has been lost. It could be argued that (3.3) 
should be regarded as its own deformation as it admits a dual (classical, and quantum) description: 
It is a solitonic system (2.1) with 

A = Hx lk' b = Hy ^x 

or a dispersionless limit (2.2) with A = XH X , B = XH y . 

One attempt to find a dispersive analogue of (3.2) would be to use the centrally extended Virasoro 
algebra in place of diff(X). Recall that such a procedure has been used to produce dispersive systems 
from dispersionless systems in a different context. Namely [23], one can view the periodic Monge 
equation ut = uu x as the equation for affinely parametrised geodesies with respect to the right- 
invariant metric on Diff(5 1 ) constructed from the L 2 inner product on the Lie algebra. Going to the 
central extension, one finds that affinely parametrised geodesies on the Virasoro-Bott group correspond 
to solutions of the KdV equation. (For a recent review of such constructions, see [21].) In the current 
situation, we view a general element of the extended algebra as a pair 

(/,«) := f(x) - iac 

where a£l does not depend on x, and c is a constant. Assuming that x is a periodic variable the 
modified commutation relation are 

%c f 

[fd x ,gd x ] c =< f,g > d x + — / (f xxx g - fg xxx )dx, 

and we see that the central term is a function of (t, y) only. Applying this procedure to (3.3) with 

H(x,y,t) = Y,h k (y,t)L k , (4.1) 

k 

where L k are generators of the centrally extended Virasoro algebra satisfying 

[L\ L m ] = (k - m)L k+m + ^k(k 2 - l)S k ,-m (4.2) 

would modify only one equation in the infinite chain of PDEs for the functions h}~. 

In the remaining part of this paper we shall present a construction 2 which leads to non-trivial 
dispersive analogues of (3.2), or its equivalent form (3.3), and (3.7). 



2 An alternative approach which we have not explored would be to consider a quantum deformation of the Virasoro 
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4.1 Deforming Lie algebra homomorphisms 

To find a non-trivial deformation one would wish to deform the Lie algebra of vector fields on S, but 
this algebra is known not to admit any non-trivial deformations [17]. 

We shall choose a different route [24], and deform the standard homomorphism between diff(S) 
and the Poisson algebra on T*£, the point being that the homomorphisms between Lie algebras can 
admit non-trivial deformations even if one of the algebras is rigid. A deformation of (3.3) is achieved 
in two steps, each introducing a parameter. In the first step we shall deform the embedding of diff(S) 
into the Lie algebra of volume-preserving vector fields on T*E. This introduces the first parameter 
pi. The second step is a deformation quantisation of the first one: The Poisson algebra on T*S is the 
quasi-classical limit of the Lie algebra of pseudo-differential operators on S, so (working at the level 
of symbols) one quantises the deformed homomorphism by using the deformed associative product of 
symbols of pseudo-differential operators rather than a pointwise commutative product of functions. 
This introduces the second parameter e. In what follows we shall be interested in the polynomial 
deformations rather than the formal ones. 

The standard embedding ir : vect(S) — > C°°(T*S) is given by contracting a vector field Xf = 
f(x)d x with a canonical one-form B on T*S. In our case T*E = K x S and the Lie algebras algebras 
C°°(S 1 ) (with the Wronskian bracket) and vect(S' 1 ) (with the Lie bracket) are isomorphic so we can 
regard ir as defined on C°°(£). 

If A is a local coordinate on the fibres of T*S, and 6 = Adx, the map ir is explicitly given by 

(7r(/))(A,x) :=A/(x). 

It is a Lie algebra homomorphism as 

Mf),*(g)} = << f,9>)- 

Given /ieM define [24] 

M/))(A,x) = (ir(f))(\,x + t x/\) = \f(x + t x/\) 

= \(f(x) + f'(x)^- + ±f"(x)^) 2 + ...). 

Note that {7T M (/), n^g)} = ^(< f,g >), so that ir^ is also a Lie algebra homomorphism between 
diff(S) and sdiff(T*£). 

algebra as in [28] and its free boson realisation 

[T m ,T n ] = — ^ /;(r n _ ; r m+ ; — T m _ ; T n+; ) — -(s™ — S m )<5 m +n,0, (4.3) 

1=1 1 ~ S 

where s = qt~ x , and coefficients /; are given by 

/(*) = £> =^{22n rr^ z ■ 

Z — v n— 1 / 

The ordinary Virasoro algebra (4.2) is recovered as q — * 1. Applying (4.3) to (4.1) would lead to a q-deformed analog of 
(3.3). 
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The next step is motivated by the canonical quantisation A — ► d/dx. For any functions F,G € 
C°°(T*S) which are also allowed to depend on a parameter \x define the Kupershmidt-Manin product 

F*G = Y- — — (44) 
^ k\ d\ k dx k ( ] 

k=0 

(this is equivalent, under an e-valued change of variable to the Moyal product) and set 

{F,G} £ = -^(F*G-G*F). (4.5) 

The Poisson bracket is recovered in the limiting procedure 

hm{F,G} e = {F,G}, 

e— >0 

but the deformed bracket is equal to the Poisson bracket for all e if F, G are linear in A. This is why 
the first deformation parameter (i is needed. 

We are now ready to propose the dispersive analog of the dispersionless equation (3.3) and (3.7). 

1. Let H(\,x,y,t; fj,,e) = ir^H) take values in an algebra of formal power series in e with an 
associative product defined by (4.4). The deformed analogue of equations (3.3) is: 

H x t - H yy — {H x , H y } 6 = 0. (4.6) 

2. Let Ni(\, x l ; /j,, e) = ^(iV;) take values in an algebra of formal power series in e with an associa- 
tive product defined by (4.4). The deformed analogue of equation (3.7) is: 

d x Ni + e ijk djN k - ^e ijk {Nj ,N k } e = 0. (4.7) 

Given a solution H of (4.6) such that (/xd„ — Xd\)(X~ 1 H) = 0, and A" 1 ^ is smooth in (///A) we 
can construct H(x, y, t) satisfying (3.3) by taking any of the two limits \x — ► 0, e — > and similar 
remarks hold for equation (4.7). Conversely, formal powers series (in the deformation parameters) 
solution may be constructed from a solution to the original, undeformed, equation in an analogous 
manner to the way developed in [29]. The extent to which such formal series converge in a suitable 
space of functions is as yet, however, unclear. 

These deformed equations formally retain their integrability; the various manipulations hold at 
the level of the Lax pair as well as at the level of the equations themselves. However, as remarked 
in the introduction, a direct twistor theory correspondence for these equations is lacking, though one 
should be able to adopt the methods developed by Takasaki [33] and FormanskL-Przanowski [9, 10] to 
study the geometry of the corresponding Riemann-Hilbert problem (in some suitable Moyal algebra 
valued loop group). 
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